The ponderomotive force due to the intrinsic spin in extended fluid and kinetic models 
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In this paper we calculate the contribution to the ponderomotive force in a plasma from the electron spin 
using a recently developed model. The spin-fluid model used in the present paper contains spin-velocity cor- 
relations, in contrast to previous models used for the same purpose. Is its then found that previous terms for 
the spin-ponderomotive force are recovered, but also that additional terms appear. Furthermore, the results due 
to the spin- velocity correlations are confirmed using the spin-kinetic theory. The significance of our results is 
discussed. 
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I. INTRODUCTION 

The ponderomotive force plays a crucial role in the non- 
linear dynamics of plasmas. Phenomena induced by thepon- 
deromotive force include, e.g. wakefield generation (ll |2|], 
soli ton formation, self-focusing and wave collapse [3|]. The 
classical expression for the ponderomotive force in a magne- 
tized plasma first derived by Karpman et al |4] was recently 
generalized to account for quantum mechanical effects, in par- 
ticular due to the electron spin ISfl. The physics of quantum 
plasmas (for recent reviews, see e.g. 0]) has recently re- 
ceived much interest due to applications in, e.g. quantum 
wells |8], plasmonics |9] and spintronics fTo*]. In particu- 
lar, the dynamical effects due to the intrinsi c sp in of the elec- 
tron has been investigated using both fluid lllll- ll3|] and semi- 
classical kinetic approaches |[l4l - jl8ll . In Ref. \y\ the pon- 
deromotive force related to the magnetic dipole moment of 
the electrons (due to the spin) was shown to induce a spin- 
polarized plasma, i.e. the spin-up and spin-down states of the 
electrons are separated. This effect was shown to be pro- 
nounced also for a plasma of relatively modest density ||5|], 
and the expressions has recently been applied to astrophysi- 
cal plasmas lfl9l l2^ . Furthermore, the spin-polarization, in 
turn, was shown to induce cubically nonlinear terms that may 
influence the high-frequency dynamics L5J. 

When the spin contribution of electrons to the ponderomo- 
tive force was calculated in Ref. f^, a relative simple fluid 
model was used. The model included the magnetic dipole 
force, spin-precession as well as the magnetization current 
due to the spin, but not the spin- velocity correlations. How- 
ever, in a recent work ||2Tll a fluid model of spin that includes 
spin-velocity correlations was shown to capture the spin ef- 
fects of kinetic theory much more accurately, although the 
comparison between fluid and kinetic theory in Ref. IuHi was 
limited to linear phenomena. 

In the present paper, we aim to improve the calculation 
made in Ref. lUt] to capture the effects of the spin- velocity 
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tensor using the newly developed four-equation hierarchy 112111 
for spin-fluid dynamics. Furthermore, in order to validate this 
model, we compare our results with the full kinetic treatment 
from the equations derived in Ref. |22]. The set-up is that of 
a weakly nonlinear electromagnetic pulse propagating paral- 
lelly to an external magnetic field. It is then confirmed that 
a spin-ponderomotive term of the same kind as in previous 
calculations exists (together with the classical ponderomotive 
force), but it is also found that the spin-velocity correlations 
induce an additional term. This new term, due to the spin- 
velocity correlations, gives rise to a force in the same di- 
rections for all particles, regardless of the spin-states (up or 
down) relative to the external magnetic field. This is in con- 
trast to the previously found term |5] which acted in oppo- 
site directions for spin-up and down states. Furthermore, the 
new term has higher order resonances at the spin precession 
frequency, as compared to the previous term |5]. The com- 
parisons with kinetic theory give a perfect agreement in the 
low-temperature limit. This comparison also includes a cal- 
culation of the low-frequency magnetization induced by the 
spin-ponderomotive force. The main purpose of the present 
paper is to validate the newly developed model against kinetic 
theory also in the regime of nonlinear perturbations. How- 
ever, we also point out that the newly found contribution to 
the spin-ponderomotive force is likely to be important for the 
nonlinear evolution in quantum plasmas. 



n. FLUID DESCRIPTION 



In Ref. 11211] (see also Ref. 112311 ) a fluid moment hierarchy 
was derived from a quantum kinetic theory with spin, extend- 
ing previous hydrodynamical models for particles with spin- 
1/2. The time evolution equations are given by 
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[d, + Vkdk)Pij = -PikdkU j - PjkdkUi - PijdkUk 

H ^imnPjm^n H ^imnPim^n 

m m 
+ ^^ikdjBk + ^^jkd^Bk - dkQijk (4) 
((9f + Vkdk)T.ij = - I.ijdkVk - ^ikdkVj - PjkdkSi 
+ —£ikl^ikBl + -^EiklLkiBt 

m n 
+ ^djB, - J^SiSkdjBk, (5) 

where n is the particle number density, v is the fluid velocity, q 
is the particle charge, m is the mass, P,y is the pressure tensor. 
Si is the spin density which is defined to have length h/2 for a 
coherent spin state, and Qijk is the heat flux tensor. The mag- 
netic moment of the particle is given hy jj, — gqh/Am, where 
g is the g-factor of the particle in question {g = 2.0023 for 
electrons). Summation over repeated indices i,j,k - ■ ■ — 1,2,3 
is understood and we have used the notation dt ~ d/dt and 
di = d / dxi and Eijk is the antisymmetric Levi-Civita tensor. 
The tensor E,, is the spin-velocity tensor which correlates spin 
and velocity and captures, in a macroscopic description, how 
different forces acts on different parts of the spin probability 
distribution. In Ref. [21] inclusion of was shown to re- 
produce some of the more subtle effects predicted by kinetic 
theory in the hnear regime. Note that quantum mechanical 
effects associated with particle dispersion (e.g. the Bohm po- 
tential) could have been included in the above model. How- 
ever, since such effects will not influence the ponderomotive 
force of transverse waves propagating along a magnetic field 
which is our problem of consideration, we do not include 
such terms here. 

Next, we consider circularly polarized electromagnetic 
waves propagating parallel to an external magnetic field, Bq — 
Bqz, and use the following ansatz 



2 

B = l 



E(z, f + E* (z, f 

B(z, f + B* (z, f 



(6) 



where the amplitudes are assumed to vary much slower than 
the exponential factors, and the star denotes complex conju- 
gates. Since the basic wave modes propagating parallel to Bq 
are either left- or right-circularly polarized, we have E,B 
X ± /y. Furthermore, all perturbations are small, such that 
weakly nonlinear perturbation theory is applicable. The equi- 
librium density and spin density will be denoted by no and Sq, 
respectively. 



A. Spin-Ponderomotive Force 

When thermal effects are small, the classical (superscipt 
cl) low-frequency (subscript (If)) ponderomotive force corre- 
sponding to the ansatz (|6|l (with transverse fields) is 
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where (Oc = qBa jm is the electron cyclotron frequency 1124 
(We will here limit ourselves to electrons, although all results 
in Sec. II can straightforwardly be generalized to any particle 
species fulfilling Eqs. ([TJ-llSll). Equation (|7]) was first derived 
by Karpman et al |4], and has been verified by many subse- 
quent authors. On the other hand, the spin-dependent part of 
the ponderomotive force was recently calculated by Brodin 
et al 1 5], but starting from a model without the spin-velocity 
tensor. Now, we follow their outline to calculate the spin con- 
tribution to the ponderomotive force, but taking also the parts 
originating from into account. 

To find the spin dependent part of the ponderomotive force, 
we first linearize the Eq. (|5]l and note that the only components 
of the spin- velocity vector that have a driving term are E13 and 
E23 and thus we define E± = E13 ± /E23, B± — B^i iBy and 
S± — Sx ± iSy. Neglecting the slow derivatives, i.e. derivatives 
acting on the amplitudes, we get the linear solution to Eq. (|5]) 
given by 



E+ = - 



hllriQk 



B4 



(8) 



where cOcg = IfiBo/h is the spin-precession frequency (which 
is close to the cyclotron frequency for electrons, as \g\ — 
2.0023). Iterating by plugging this back into Eq. (|5]l we find 
the correction due to the slow derivatives as 
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This last step could also be done by simply making the substi- 
tutions 0) — > O + idt and k ^ k — id, in Eq. (O and expanding 
to lowest order in the slow derivatives. From Faraday's law 
the electric and magnetic fields are similarly related by 



, k 1 dE± k dE± 

B± = ±i-E± ± — -i ± _— 

CO CO dz CO^ dt 



(10) 



The lowest order approximation, B± = ±i{k/co)E±, can be 
used to switch between magnetic and electric fields in, e.g. 
the right-hand side of Eq. O. Repeating the above steps for 
Eq. (|3]l and using Eq. (|9]l we can express the spin variable as 
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The spin-ponderomotive force density is then given by the last 
term of Eq. ([3]) and can be written as 



-Sld,B±). 



(12) 



It should be noted that the pressure tensor in Eq. dH actually 
has quadratically nonlinear source terms proportional to the 
magnetic moment which, in principle, could contribute to a 
low-frequency spin force, in addition to that originating di- 
rectly from the magnetic dipole force. However, following 
the same calculation scheme as outlined above, it can be ver- 
ified that these terms do not contribute to the leading order in 
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the slow derivative expansion. Thus, using Eq. ( fTTT ) to the first 
order in slow derivatives we obtain 
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B. Low-frequency magnetization 

The slowly modulated wave also gives rise to a nonlinearly 
induced low-frequency magnetization. To further investigate 
the fluid model we next calculate the low frequency magneti- 
zation defined as 



Mif = ^(noSif + nifSo) 



(14) 



Comparing this force with the result of Ref. we see that 
the first terms (proportional to ^o) corresponds exactly to 
their result [after a correction of factor 2 in Eq. (13) of Ref. 
pt] has been made]. The second term comes in due to the 
inclusion of the spin-velocity tensor (which was neglected 
in Ref. [5]). It has an extra factor Tik^ /2m{(0^ (Ocg) com- 
pared to the first term and could in a sense be viewed as 
a higher-order quantum correction. This term might, how- 
ever, dominate over the first term in some cases where the 
zeroth order magnetization is small. For thermal equilibrium 
we have Sq = {ti/2)\.imh{}iBQ/kBT) where tanh(^BoABr) 
~ jJ-Bo/ksT for a moderate magnetic field strength, and thus 
both the terms of the spin-ponderomotive force are quadratic 
in h. In this case, the ratio of the first to second term scales 
as ^ m(Ocg (ft) — (Ocg) /ksTk^, which could be both smaller or 
larger than unity. On the other hand, for a higher magnetic 
field strength for which tanh (/iBq I^bT) ^ 1 , the same ratio is 
scaled as ^ m (ft) — (Ocg) /hk^ , which can also be either smaller 
or larger than the unity. Thus, depending on the parameter 
regimes we consider for a specific problem, both the classical 
and the spin-ponderomotive force may either be comparable 
or even dominate over each other lfT9ll20ll . As for example, for 
a wave frequency close to the electron-cyclotron frequency, 
the spin contribution to the pondero motive force can, indeed, 
dominate over the classical one when hk^/m(0 > 1 filillH. 

A slightly different perspective was taken in, e.g. Ref. lit]. 
There a two-fluid model for electrons was used, where spin-up 
and down states were treated as different species, and thus the 
two species have Sq — ±{h/2). As a consequence, the force 
on each separate species [due to the first term of Eq. (fT3] )1 is 
typically stronger [since the diminishing factor hBq /ksT dis- 
appears from the first terms, the ratio between the two types of 
terms are now hk^ /2m (ft) ^ (Ocg)], but on the other hand, the 
forces on the two spin-populations act in opposite directions 
due to the terms proportional to Sq. Since the physics is differ- 
ent depending on whether the terms induce spin-polarization 
or not, both types of terms may be needed for an accurate 
description even in the cases when one of the terms is much 
larger than the other As argued in Ref. ll2lll . the two-fluid 
model captures some of the effects of a kinetic theory, but it 
is less needed when the spin-velocity correlations of the 
tensor is included in the model. We will return to this issue 
of one-fluid versus two-fluid models in a little more detail be- 
low. As a final note of this subsection we point out that for 
slowly modulated waves as considered here we may in many 
cases use the approximate relation (9, k, —Vgd, (where Vg is the 
group velocity) to simplify Eqs. Q and (fl3] l, although more 
accurate expressions might be needed in case the terms from 
the spatial and temporal derivative tend to cancel. 



We are particularly interested in the case where the low- 
frequency magnetization is induced in the absence of prior 
magnetization, and hence we take So = for simplicity. As 
Sq — {h/2) tanh {iJ.Bo/kBT) in a thermal equilibrium back- 
ground, this may be a useful approximation also from a prac- 
tical perspective, provided jiBo/kBT ^ 1. 

Applying Eq. (O, considering the low-frequency time scale 
and keeping terms up to the second order in amplitude, we 
obtain 



d,S{, = T^^{iBiS*^ 



-c.c. 
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where c.c. denotes the complex conjugate. The first term 
can be calculated from Eq. (TT[ . and we note that only terms 
involving slow derivatives survive when adding the complex 
conjugates. As for the second term, only the Ejj components 
need to be calculated. Thus, out of these three, the only com- 
ponent with a nonzero driving term is 



dt^'^±i-^i:±B*^+c.c. 

In 



(16) 



where +c.c. means adding the complex conjugate of the right- 
hand side. Combining the above formulas with the linear re- 
sults from the previous section, Eqs. (|9]l and ( fTTI ) we obtain 
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C. Magnetization in spin two-fluid models 

In Ref. f3\ a similar problem was studied, where the first 
term of Eq. ([T3T l induced opposite density perturbations for 
particles with initial spin-up (i.e. with5o = +ti/2) and for par- 
ticles with initial spin-down (with = —h/2)). As pointed 
out above, it is straightforward to include such a formalism 
within the present theory, we just consider two electrons flu- 
ids which are identical in all aspects, except that the unper- 
turbed spin in the z-direction is +h/2 or —Ji/2. The ex- 
pression Eq. ( [T3] ) is then the same, except that we substi- 
tute So = ±{h/2) where + (— ) stands for up (down) species, 
rather than the thermodynamic equilibrium expression = 
{h/2)tanh{jiBo/kBT) of the one-fluid theory. 

Let us now calculate the two-fluid version corresponding 
to Eq. fm . which is still evaluated when fiBo/kBT ^ 1 such 
that of the one-fluid theory is negligible. This corresponds 
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to equal densities of the initial spin-up and down- popula- 
tions in the two-fluid model. Somewhat surprisingly we ob- 
tain exact agreement with Eq. (l% using such a two-fluid the- 
ory. However, contributions to the magnetization that come 
from the spin-velocity correlations in the one-fluid model is, 
to some extent, replaced by contributions that come from den- 
sity perturbations of the two-species in the two-fluid model. 
Specifically, we note that our current model (which, as we 
recall, includes a modified ponderomotive force due to the 
spin-velocity tensor) agrees with the density perturbations of 
the spin-up and spin-down populations of the less elaborate 
model used in Ref. |5]. Furthermore, we stress that the ad- 
ditional spin-ponderomotive force term does not contribute to 
the spin-polarization. Nevertheless, the corresponding mag- 
netization found here can be much lower in certain cases, as 
compared to the model without spin-velocity correlations. In 
particular, this holds for the specific case of an unmagnetized 
plasma {Bq — 0), and a weakly dispersive driver, i.e. from 
Eq. ( [TtI i we find that for d, = — v^^., Bq — O and zero disper- 
sion (i.e. Vg = co/k), the right-hand side of Eq. ( [TtI i vanishes. 
This may seem surprising, as intuitively one would expect that 
a net density difference between the up- and down-species 
(as found in a two-fluid model both with and without spin- 
velocity correlations) implies a net magnetization. However, 
the reason of why this does not give rise to a finite magneti- 
zation for the given assumptions specified above, is that the 
density difference is compensated by a rotation of the spin di- 
rection, which also is of second order in the transverse wave 
field amplitudes. 

In a two-fluid model [where the unperturbed spin of the 
species are Sq = ±{h/2)z] the magnetization for each species 
is then given by Eq. (fl4] i and the total low-frequency magne- 
tization is obtained as 



2jU 



h 



no('Szifr - Sz\fi) + (nift - n\fi)So] 



(18) 



with So — h/2. For the assumptions specified above where the 
induced magnetization vanish, we accordingly have (^^if^ — 
= -("ift -«ifj,)^V2no- 



III. KINETIC APPROACH 

In order to further investigate the fluid model of Eqs. 
(|5]l, we now consider the same situations as above but using 
a kinetic model. The evolution of the scalar distribution for a 
spin particle is, in the long scalelength limit, given by t22il 

^r/ + v-V,/+-^(E + vxB)-V,./ 
m 

+ ^V,- (B -s + B • V,) ■ V„/+ ^(s X B) • V,/ = 0. (19) 
m n 

Here, the distribution function / = /(x, v,s,f) contains a spin 
variable s, which is defined to lie on the unit sphere. This 
can be directly related to the fluid model, with the macro- 
scopic spin vector given by S(x,f) = {h/2) ^ dQ.3sf{x,\,s,t), 



where the integration element is dCl ~ cfivcfis and the two- 
dimensional spin integration is carried out over the unit 
sphere. For a more detailed description we refer to Ref. ll22ll . 
However, we point out that the spin part of the distribution 
function is more general than a semi-classical treatment, such 
as that made in, e.g. Refs. lfl4llT5ll . where the evolution equa- 
tion for the distribution function has the same structure as for 
a classical magnetic dipole moment. 

In order to calculate the weakly nonlinear low-frequency 
response to an incoming transverse wave packet we make the 
ansatz 

/(x,v,s,f) =/o(v^ 0,) +/if(z,f ,v,0,) 

(20) 

where /o is the background distribution, /if is a low-frequency 
part due to quadratic nonlinearities and f\ is a slowly modu- 
lated high-frequency wave. The background distribution will 
be taken to be of the form 



/o 



no 



27r3/2v3 



(21) 



where «o is the eq uilibrium density, the thermal velocity vj is 
defined as vj = ^/2Jq}T/m. This is the generalization of ther- 
mal equilibrium for the scalar distribution function, applicable 
for low or moderate densities, which contains the quantum 
effects on the angular spin distribution. The fully quantum 
mechanical background distribution (applicable in the regime 
of high densities and/or strong magnetic fields) is presented 
in Ref. [22]. The spin part of the distribution function above 
can also be written as f^{\ + cos 0s) +/;(1 — cos0), where 
/j- exp{-^BBQ/kBT) and exp{iiBBo/kBT). Thus, we 
see that the part proportional to cos 9s scales as (/| — fi)/{f-\- + 
fl), which gives the factor tanh{iiBBo/kBT) in the second 
term of Eq. (l2Tl i. It is then more clearly seen that the distri- 
bution function is separated into two parts, one for each spin 
population. This can be considered as a basis for derivation 
of a two-fluid model (2A)i, but we will not pursue this further 
here. Note also that a semi classical reasoning would suggest 
us to take the distribution function as /o exp[—E/{kBT)], 
where E — mv^/2 — jiBocos 9s, which differs from the quan- 
tum mechanical angular distribution used here. We will as- 
sume that the temperature is sufficiently low for vj- ^ co/k 
such that we can neglect thermal effects in the final result. 

The ansatz for the incoming high-frequency field in Eq. ^ 
is the same as before. The aim is then to find an equation for 
the low-frequency part of the distribution function. From such 
an equation we can then calculate the low-frequency response 
in the current density and magnetization, and compare with 
the results from the previous section. The high-frequency per- 
turbation of the distribution function is fi = (/_ ) for left- 
hand (right-hand) circularly polarized waves. The expressions 
for f± found from Eq. (fT9] l to linear order can be easily com- 
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puted from Eq. (9) in Ref. |125| 



f± 



(0 — kv^ =F Oc 2m 



E±dy^fo 



(22) 



2m 



kB± (sin Osdv-fo + cos Osde^dv-fo) ± —B±deJo 



Next, allowing for slow modulations and solving the equation 
to first order in d^/k, d, / co, we note that the zero-order so- 
lution applies after making the substitution O) — > 0) + id, and 
k—7'k — id, in Eq. (|2TI) . and then expanding to first order in the 
slow derivatives. Inserting the ansatz above into the evolution 
equation and consider the slow-time scale and keeping only 
up to quadratic nonlinearities we obtain the equation 



{d, + v,d^) fif + ^E,ifdyJo 

(E + V X B) + i^- (/fc + d,) (s • B + B • V,) zl • V,./i* 
Am ^ ' 4m ^ ' \ ' 



^sxB-Vjf+c.c. 



(23) 



Here we have also added a low-frequency electric field in the 
z-direction, fi^if, which has /if as source. Equations (l22l i and 
(l23T l now constitute a basis for calculating the nonlinear re- 
sponse in the current density and magnetization. 



A. Ponderomotive Force 

Inserting the first order solution /i in Eq. (|23] |. multiply- 
ing by qv- and integrating over d£L = d^vcfis we can derive 
an equation for the current density J,. We will neglect Lan- 
dau damping associated with the particle resonances. Further- 
more, since we have assumed a low-temperature (vj- ^ (£> jk) 
we may expand the denominators in f\ to the lowest order in 
v^. All the integrals can then be evaluated. 

As an intermediate step, after integrating Eq. ( l23T l we obtain 

2 

dtJ^i + ?<5' / dQ.v}f]f — - — -E,\f = nonlinear terms. (24) 
"J ' m 

The second term on the left-hand side of Eq. ( |24] | is a thermal 
correction to the low-temperature limit, which we will neglect. 
Using Ampere's law for the first term on the left-hand side, 
we can then obtain a closed equation for the time derivative 
of the low-frequency electric field in terms of the incoming 
wave. After some algebra we find 



[d^ + coj]E,ii^C±+D±, 



deromotive source terms on the right-hand side are given by 

cOck 
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Note that the term C± is due to the classical ponderomotive 
force stemming from the magnetic part of the Lorentz force, 
and D± is due to the spin effects. Taking into account the 
normalization factor no /q (which appears when an evolution 
equation for E^if is derived) we have perfect agreement with 
the fluid theory in the previous section. 



B. Magnetization 

We can also obtain the low-frequency response in the mag- 
netization. This is done by multiplying Eq. ( |23l ) by 3/1^5- and 
integrating over the velocity and the spin. Thus, we obtain 



dtM^\f + 3nd, / (iHsj-v^/if = nonlinear terms. 



where cOp — noq^ / meo is the plasma frequency and the pon- 



(27) 



In order to be able to evaluate the second term, we take the 
time derivative of Eq. ((|27]l) and once again use Eq. ( |23] l to 
evaluate the resulting term 3/iB(9j/ dQ.s^v,dtf\f. We still ne- 
glect all the thermal contributions including the particle res- 
onances. To simplify the problem further, and compare with 
the fluid result, we also neglect the zeroth order magnetiza- 
tion which corresponds to neglecting the factor proportional 
to = {h / 2)x.w\\{iiBQ / kaT) in Eq. dSTT i. Also, to obtain an 
evolution equation for Mj-if, one takes the time derivative of 
Eq. ( |27] i, substitute ( |23] l for (9,/if, apply Eq. ( l22l i and keep 
only the low-frequency source terms (those proportional to 
or liiip), and carry out the integrals over dQ., omitting 
thermal corrections. After lengthy calculations and evaluat- 
ing the corresponding integrals we eventually obtain an exact 
agreement with Eq. ( fTTl ). This is an important verification 
that the truncation of the moment hierarchy used to close the 
fluid equations is valid in the low-temperature limit, also when 
nonlinearities are present. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we have calculated the ponderomotive force 
due to the electron spin property in both a recently derived 
fluid model lf2lll and in the kinetic model (see Eq. (83) in Ref. 
['22']) that is the basis for the fluid theory. The kinetic result 
has been evaluated in the low-temperature limit, in which case 
we obtain a complete agreement with the fluid result. The 
fluid theory considered here extends a simpler fluid model 
(e.g. Ref. jsll) by including the spin-velocity correlations. As 
a result, the spin part of the ponderomotive force found here 
has a contribution that is additional to that recently derived 



6 



in Ref. (3}. The previously derived spin force contained a 
term proportional to the unperturbed spin So (the first term in 
(fTsl i) but not the second term in Eq. ( fTSI l that is independent of 
So- Nevertheless, Ref. |5] found that the spin-ponderomotive 
force could be important even when the unperturbed magne- 
tization (and hence the net value of So) is zero, such as in an 
unmagnetized plasma. The reason was that a two-fluid model 
of electrons was used, where the spin-up and spin-down states 
were exposed to spin-ponderomotive forces pointing in oppo- 
site directions, which induced a spin-polarized plasma. It is 
straightforward to include such an approach also within Eqs. 
([TJ-llSll. Furthermore, within a nonlinear perturbation scheme, 
the division of the electrons into spin-up- and spin-down- 
populations can be put on a firm basis in the more complete 
kinetic theory. 

However, provided the more subtle effects of spin-velocity 
correlations are taken into account, our results here, as well 
as those in Ref. ||2lll , suggest that a two-fluid spin model 
(which may capture certain effects of the microscopic spread 
in the spin-probability distribution) is not needed. This results 
from the fact that the spin-velocity tensor seems to capture 
the physics of spin-polarization, as well as additional effects 
of the spin. On the other hand, it is still too early to exclude 



the possibility that treating spin-up and -down populations as 
different species captures new physics compared to a one-fluid 
model. In particular, a difference between the one- and two- 
fluid models may reveal itself when cubically nonlinear cal- 
culations are carried out. 

One of the main conclusions of this paper is that Eqs. ([TJ- 
(|5]l agree completely with kinetic theory when thermal effects 
are negligible, whereas the same would not be true when the 
spin-velocity tensor is omitted. This comparison holds not 
only for the spin-ponderomotive force, but also for the in- 
duced low-frequency magnetization, that has been calculated 
both from the fluid and the kinetic theory. It should be noted 
that the previous comparison between the current spin fluid 
model and the spin kinetic theory was limited to linearized 
theory. The previous expressions for the ponderomotive force 
were shown to be significant in plasmas of relatively modest 
density and also when the plasma was unmagnetized, due to 
the induced spin-polarization f?]. The additional term in the 
expression for the spin ponderomotive force found here may 
be at least as important, particularly when the wave frequency 
is close to (Ocg, due to the higher-order of the resonance in 
this new term. A detailed evaluation of the effect of the spin- 
ponderomotive expression found here is a project for future 
work. 
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